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Abstract: In this paper we established the condition for a curve to satisfy stochas- 
tic fractional HP (Hamilton-Pontryagin) equations. These equations are described using 
Ito integral. We have also considered the case of stochastic fractional Hamiltonian equa- 
tions, for a hyperregular Lagrange function. From the stochastic fractional Hamiltonian 
equations, Langevin fractional equations were found and numerical simulations were done. 
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1 Introduction 



in 
o 

g . J.M. Bismut was the first one that introduced concepts of stochastic geometric me- 

chanics, in his work from 1981, when he defined the notion of "stochastic Hamiltonian 
system" . He showed that that the stochastic flow of a certain randomly perturbed Hamil- 
^ I tonian systems on fiat spaces extremizes a stochastic action, and using this property, 

■ he proved symplecticity and Noether theorem for stochastic Hamiltonian systems. Since 

then, there has been a need in finding out tools and algorithms for the study of this kind 
of systems with uncertainty. Bismut's work was continued by Lazaro-Cami and Ortega 
(|17j. [T8]). in the sense that his work was generalized to manifolds, stochastic Hamilto- 
nian systems on manifolds extremize a stochastic action on the space of manifold valued 
semimartingales, the reduction of stochastic Hamiltonian system on cotangent bundle of 
a Lie group, a counter example for the converse of Bismut's original theorem. 

Very important in many science fields is fractional calculus: fractional derivatives, 
fractional integrals, of any real or complex order. Fractional calculus is used when frac- 
tional integration is needed. It is used for studying simple dynamical systems, but it 
also describes complex physical systems. For example, applications of the fractional cal- 
culus can be found in chaotic dynamics, control theory, stochastic modelling, but also 
in finance, hydrology, biophysics, physics, astrophysics, cosmology and so on ([H], [S], 
[in]). But some other fields have just started to study problems from fractional point of 
view. In great fashion is the study of fractional problems of the calculus of variations 
and Euler-Lagrange type equations. There were found Euler-Lagrange equations with 
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fractional derivatives, and then Klimek found Euler-Lagrange equations, but with sym- 
metric fractional derivatives [16]. Most famous fractional integral are Riemann-Liouville, 
Caputo, Grunwald-Letnikov and most frequently used is Riemann-Liouville fractional in- 
tegral. The study of Euler-Lagrange fractional equations was continued by Agrawal [2] 
that described these equations using the left, respectively right fractional derivatives in 
the Riemann-Liouville sense. This fractional calculus has some great problems, such as 
presence of non-local fractional differential operators, or the adjoint fractional operator 
that describes the dynamics is not the negative of itself, or mathematical calculus may 
be very hard because of the complicated Leibniz rule, or the absence of chain rule, and so 
on. After O.P. Agrawal's formulation of Euler-Lagrange fractional equations, Baleanu 
and Avkar [4J used them in formulating problems with Lagrangians linear in velocities. 
Standard multi-variable variational calculus has also some limitations. But in ^23] C. 
Udri§te and D. Opri§ showed that these limitations can by broken using the multi- linear 
control theory. 

For fractional stochastic integrals of the form 

Xt = Xo+ [ F{t,s,X,)ds+ [ Ga{Hut,s,Xs)dW:, (1) 
Jo Jo 

the existence and uniqueness of its solution was discussed by Pardoux and Protter in their 
work [2D]. 

In , it was proved the existence, uniqueness, and continuity of a fractional stochastic 
equation of the form 

Xt = Xo + I^F{t, s, X,) + WfG^iHf, t, s, X,), (2) 
where < /3 < 1, It F{t^ s, Xs). The fractional integral of F{t, s, Xs) is defined by 

I^Fit sX)-^t ^:^^i^ds 

and fractional Wiener process WfGaiHf, t, s, Xg) of GaiHf, t, s, Xg) is 



One application of fractional stochastic equation of the form ^ is in finance. Frac- 
tional Black-Scholes market is described in terms of the bank account and a stock. The 
price at a time t is given by the following formula 

At = exp(^ J r(s)ds^, (3) 

where r(s) > 0, s G [0,t], represents the interest rate. The price can be expressed using 
a fractional Volterra-type equation: 

where cr > are continuous functions on [0, T]. 
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In this paper, we restrict our attention to stochastic fractional Hamiltonian systems 
characterized by Wiener processes and assume that the space of admissible curves in 
configuration space is of class C^. Random effects appear in the balance of momentum 
equations, as white noise, that is why we may consider randomly perturbed mechanical 
systems. It should be mentioned that the ideas in this paper can be readily extended to 
stochastic Hamiltonian systems [19] driven by more general semimartingales, but for the 
sake of clarity we restrict to Wiener processes. Within this context, the results of the 
paper are as follows: 

1. The paper proves almost surely that a curve satisfies stochastic fractional HP 
equations if and only if it extremizes a stochastic action. This theorem is the main result 
of the paper; 

2. Fractional HP equations are described using fractional Riemann-Liouville integral 
and fractional Ito integral; 

3. Langevian type stochastic fractional equations are obtained, in the case of a hyper- 
regular Lagrange function. 

The paper is organized as follows. In Section 2 we present some sufficient conditions 
for existence, uniqueness and almost sure differentiability of stochastic flows on manifolds. 
In Section 3, we extend the fractional Hamilton-Pontryagin (HP) principle to the stochas- 
tic setting to prove that a class of mechanical systems with multiplicative noise appearing 
as forces and torques possess a variational structure. In Section 4, for a hyperregular 
Lagrange function, we get the stochastic fractional Hamiltonian equations that lead to 
Langevin fractional equations. For a fractional Lagrange function, defined on M^, the 
corresponding fractional Langevin equations are simulated. The mechanical system could 
evolve on a nonlinear configuration space and involve holonomic constraints or noncon- 
servative effects in the drift. The fractional Hamiltonian and the Lagrangian description 
are joined together to get the fractional HP system. 

2 Stochastic flows on manifolds 

Some standard results on flows of SDK on manifolds are reviewed here to the reader's 
convenience. For more detailed exposition, the reader is referred to the textbooks such 
as [12] or [15]. This section parallels the treatment of deterministic flows on manifolds 
found in [T]. 

Let M be a manifold modelled on a Banch space E. Recall that a vector fleld on the 
manifold M is a section of the tangent bundle TM on M. The set of all C'^ vector flelds 
on M is denoted by X^(M). 

A stochastic dynamical system consists of a base flow on the probability space which 
propagates the noise, and a stochastic flow on M which depends on the noise. 

A stochastic dynamical system consists of a base flow on the probability space {Q, 3", P) 
and a stochastic flow on a manifold M. The base flow is a P— preserving map 9 : M.xfl ^ Q 
which satisfles: 

1. 9q = icLq : VL ^ VL is the identity on 

2. for all s, t G M, the group property,^^ 06^ = 6t+s- 

Given times < r < s < t, the "stochastic flow" on M is a map ift^^ '■ ^ x M ^ M 
such that 
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1. for almost all a; G fi, the map {s,t,uj,x) ^ ipt,s{^)x is continuous in s,t and x; 

2. (y9s^s(a;) = idM : M — > M is the identity map on M, for all s G M; 

3. if satisfies the co cycle property 

This paper is concern with stochastic dynamical systems that come from stochastic 
laws of motion, i.e. ones whose stochastic flows define solutions of stochastic differential 
equations. Consider a manifold M, modelled on a Banach space E and vector fields 
Xo,Xa G X'^(M), a = l,...,m. Let tu), 3",-), a = l,...,m, be independent Wiener 

processes for < t < T. In terms of these objects, the Stratonovich stochastic differential 
equations, that the paper considers, takes the form: 

dx = Xo{x)dt + Xa{x) o dW, x{0) = xq. (5) 

A Stratonovich integral curve of (0) is a C"— map, c{-,uj) : [0,T] M which satisfies 

c{t,io) = xo+ / Xo{c{s,u))ds + / Xa{c{s,uj)) o dW'^{s,u)ds, 
Jo Jo 

for all t G [0,T]. 

Let c be a Stratonovich integral curve of ([5]). Pathwise uniqueness of c means that if 
c : J — M is also a solution of ([5]) on the same filtered probability space, with the same 
Brownian motion and initial random variables, then, 

P(c(t, u) = c{t, u), Vt G [0, T]) = 1. 

For the rest of the paper the explicit dependence of stochastic maps on the point a; G f2 
will usually be suppressed. With these definitions, one can state the following key, but 
standard theorem ([H], [12], [15]). 

Theorem 1 (Existence, uniqueness and smoothness) 

Let M be a manifold with the model space E. Suppose that XQ,Xa G X^{M), a = 1, ...,m 
and k >1 are uniformly Lipschitz and measurable with respect to x & M. Let I = [0,r]. 
Then the following statements hold. 

1. For each u G M, there is almost surely a C'^ — curve, c : I ^ M, such that c(0) = u 
and c satisfies ([5]), for all t. This curve c : I ^ M is called a maximal solution; 

2. c is pathwise unique; 

3. There is almost surely a mapping F : I x M ^ M,such that the curve Cu '■ I ^ M, 
defined by Cu(t) = Ft{u),is a curve satisfying (ED, for all t & L Moreover, almost 
surely F is C*^ in u and C° in t. 

□ 



4 



Chi§ Oana, Opri§ Dumitru 



3 Stochastic fractional HP mechanics 

In this section, a fractional variational principle is introduced for a class of stochastic 
fractional Hamiltonian systems on manifolds. The stochastic fractional action is a sum of 
classical fractional action and several stochastic integrals. 

Let Q be an n— dimensional manifold, 7a : Q — >■ M, a = l,...,m, a deterministic 
function, and the Lagrangian L : TQ — R. Let (fi, 5", P) be a probability space and 
the interval [a, 6] C M. Let {W°'{t),3^t\t(z[a^, for a = l,...,m, where {W^}a=i,...,m are 
independent, real- valued Wiener processes and 3^t is the filtration generated by these 
Wiener Processes. 

The stochastic HP fractional action is defined by A"''^ : Q x C[PQ) M: 
A^^^iq, v,p, t) = rfe) JaiLiQis), ^i^W - s)""' + ivis), q{s) - v{s)){t - 



s 



(6) 

where 0<a<l, 0</?<l, T{a), T{P) are the Euler gamma functions and PQ = 
TQ © T*Q, and 

e(PQ) = {{q,v,p) e C''{[a,b],PQ)\qeC\[a,b],Q), q{a) = q^, q{b)=qk}. (7) 

The first integral in (I6l) is a Riemann integral, and the second one is an Ito integral. 

The HP path space is a smooth infinite-dimensional manifold. One can show that its 
tangent space in c = {q,v,p) G C([a, 6], gi, ^2) consists of maps w = {q,v,p,6q,6v,6p) G 
C°([a, 5], T(P(5)), such that 6qa = Sqt = 0, and q, 6q are of class C^. Let us denote by 
{q,v,p){-,e) G G{PQ) the one-parameter family of curves in C, that is differentiable with 
respect to e. Define the differential of A""'^ as 

dA'''''{6q,Sv,Sp) = ^A'''''{u;,q{s,e),v{s,e),pis,e)) , (8) 



e=0 



where 



d 



de 



d 



^lis) = -7^(lis,e) , Sq{a) = 5q{h) = 0, 5v{s)—v{s,e) , dp{s) = —p{s,e) . (9) 



£=0 de 



d 



e=o de 



e=0 



In terms of this differential, one can state the following critical point condition for the 
action A"'^. 

Theorem 2 (Stochastic fractional variational principle of HP) 

Let L : TQ ^ M. be a Lagrangian on TQ of class C^, with respect to q and v, and with 
globally Lipschitz first derivative with respect to q and v. Let 'ja '■ Q ^ ^ be of class 
and with globally Lipschitz first derivatives, for a = 1, ...,m. Then, almost surely, a curve 
c = (g, V, s) G G{PQ) satisfies the stochastic fractional equations 

dq{s) = v{s)ds, 

Ms) = {^{q{s)Ms)) + §{q{s),v{s))^ys + ^^^(t - s^-^^ o dW^is), 
Pis) = §{q{s),v{s)),se[a,b], 

(10) 

if and only if it is a critical point of the function A'^'^ : Q x Q[PQ) M, i. e. dA°''^{c) = 0. 
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Proof: The proof results by applying the method from [5] and [TO]. □ 

Observe that by the Ito-Stratonovich conversion formula, the Ito modification to the 
drift is equal to 0, and hence fllOl) can be written in the Ito form as 

dq = vds, 

dp=[^{q,v) + §{q,v)^)ds+^^^{t-sr-dW^{s), (11) 

p=^{q,v), s G [a,b]. 

In what follows, structure-preserving properties of the flow map, defined by maximal 
solution of the equations over [a, b] will be investigated. First, observe that because 
of smoothness conditions assumed in Theorem 2, a solution almost surely exists and 
it is pathwise unique on [a,b], by the result from Section 2. When 7^ is constant for, 
a = l,...,m and a = 1, the reader is referred to [21], for deterministic treatments of 
symplecticity, momentum map preservation and holonomically constrained mechanical 
systems. 

If P = a, the equation (|TT1) is given by: 
dq = vds, 

dp= {^{q,v) + ^{q,v)^ys+'-^dW^{s), (12) 

If 7a is constant, for a = 1, m, from ( fTTl) . results: 

dq = vds, 

dp= [^iq,v) + §iq,v)^)ds, (13) 
and they represent the Euler-Lagrange fractional equations (^ID], |13j). 

4 Stochastic fractional equation for the Lagrangian 
hyperregular 

Let L : TQ — M be a Lagrangian on TQ hyperregular, that means det ^ q'^^q^j ^ 7^ 0. 
From ( |T0|) results the following propositions: 

Proposition 3 (Stochastic fractional Hamiltonian equations) 
The equations f fTTj) are equivalent with the equations: 

dQ' = fds, 

(14) 

dp^={-'4 + f^p)ds + ^^(t - sf-'^dW'^is), I = 1, n, 
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where H = piq^ — L{q, v). 



□ 



Proposition 4 If L = ^QijV^v^ , where gij are the components of a metric on the manifold 
Q, equations (11 II) take the form: 



dq^ = v^ds, 



dv' = -{ri.v^v' + ^v'ys + EWgijd^^t _ s)0-^dW%s), 1,3 = 1, 



(15) 



where r*;^, are Cristofel coefficients associated to the considered metric. Equations ( |T^ 
become: 

dq^ = v^ds, 



dPi = {l^P'p' + lE7P^)ds + - sr-'^dW^is), z, J, k=l,.. 



(16) 



□ 



Equations (HM represent fractional Langevin equations. Equations ( fT6i) can be used for 
fractional motion of relativistic particles with noise. 

Proposition 5 a) IfQ = M, H{p,q) = ^p'^ + U{q) and -y^q) = cos(g), equations (IT^ are 
given by: 

dq = pds, 

(17) 

dP=[-t + T^Pjd' - - sY~'^^My{q)dW{s)- 

b) IfU{q) = cos(g), the Euler scheme for (WlM s: 



x{n + 1) = x{n) + hy{n), 

y{n + 1) = y{n) + /i(sin(x(n)) + fE^y{n)) — fpy(i — nh)^'" sin{x{n))G{n) , 



where h = f , G{n) = W{{n+l)h) -W{nh), n = 0,..., A^-1, 0<a<l, 0</3<l, 
and x{n) = q{nh), y{n) = p{nh). 

For a = 0.6, P = 0.3, t = 0.8 and h = 0.0001, with Maple 13, the orbit {n,p{nh)) is 
represented in Figure 1, and the orbit {n,p{nh,uj)) is represented in Figure 2 
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Figure 1: the orbit {n,p{nh)) Figure 2: the orbit {n,p{n,uj)) 
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In Figure 3 the orbit {q{nh) , p{nh)) is represented, and in Figure 4 it is represented 
the orbit {q{nh,Lj),p{nh,uj)). 



0.398- r'---- 



1 ;-V 



Figure 3: the orbit {q{nh) , p{nh)) Figure 4: the orbit {q{nh,uj),p{nh,uj)) 

For a = (3 = 0.6, and t = 0.8, h = 0.0001, the orbits {n,p{nh, u)) and {q{nh, uj),p{nh, uj)) 
are represented in Figure 5 and Figure 6. 
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Figure 5: the orbit {n,p{nh,u>)) Figure 6: the orbit {q{nh,uj),p{nh,uj)) 



Conclusions 



In this paper it was described stochastic fractional HP principle, using classical stochas- 
tic HP principle [S] and fractional principle ([ID], [IS])- Using a hyperregular Lagrange 
function, Langevin-type fractional equations were illustrated. We have done the numerical 
simulation for the case of a Lagrange function defined on M^. 

In the future work, we will consider other problems that deal with stochastic fractional 
HP principle. 
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